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ABSTRACT 

We present results obtained from a series of observations of the supernova remnant RX J17 13.7-3946 by 
the Suzaku satellite which cover about two-thirds of the remnant surface. Hard X-rays have been detected 
from each pointing up to ~ 40 keV. The hard X-ray spectra are described by power-law functions with pho- 
ton indices of ~ 3.0, which are larger than those in the energy region below 10 keV. The combination of the 
spatially-integrated XIS and HXD spectra clearly reveals a spectral cutoff in the X-ray spectrum which is linked 
to the maximum energy of accelerated electrons emitting synchrotron radiation. The broad-band coverage of 
Suzaku observations from 0.4 keV to 40 keV allows us to derive, for the first time, the energy spectrum of 
parent electrons in the cutoff region. The inferred cutoff energy in the spatially-integrated X-ray spectrum 
indicates that the electron acceleration in the remnant proceeds close to the Bohm-diffusion limit. We dis- 
cuss implications of the spectral and morphological properties of Suzaku data in the context of the origin of 
nonthermal emission. The Suzaku X-ray and the H.E.S.S. TeV gamma-ray data together hardly can be ex- 
plained within a pure leptonic scenario, unless we introduce an additional component of relativistic electrons 
with softer energy spectrum. Moreover, the leptonic models require very weak magnetic field which does not 
agree with the recently discovered filamentary structure and short-term variability features of the X-ray emit- 
ting region. The hadronic models with strong magnetic field provide perfect fits to the observed X-ray and TeV 
gamma-ray spectra through the synchrotron radiation of electrons and p-p interactions of protons, but require 
special arrangements of model parameters to explain the lack of thermal component of X-ray emission. For the 
morphology studies, we compare the X-ray an TeV gama-ray surface brightness maps using the Suzaku XIS 
and the H.E.S.S. data. We confirm the previously reported strong correlation between X-ray and TeV gamma- 
ray emission components. At the same time the Suzaku data reveal a deviation from the general tendency, 
namely, the X-ray emission in the western rim regions appears brighter than expected from the average X-ray 
to gamma-ray ratio. 

Subject headings: acceleration of particles — ISM: individual(RX J 17 13.7-3946) — ISM: supernova remnants 
— X-rays: ISM 



1, INTRODUCTION 

Supernova remnants (SNRs) have long been considered to 
be likely acceleration sites of cosmic -ray particles below the 
energy of the knee, ~ 10 15 eV. The energy supply to ex- 
plain the energy density of comic rays is satisfied if ~ 1- 
10% of the energy of each supernova is transferred to ac- 
celerated particles. Also, the well developed theory of dif- 
fusive shock acceleration nicely ex plains the universal power- 
law spectrum of cosmic rays (e.g. iBlandford & Eichlerl lT987l : 
Malk ov~& Drurvll2001l) . Although synchrotron emission de- 
tected in the radio band supports this idea observationally, no 
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evidence of acceleration to TeV energy had been observed un- 
til recently. During the last decade, such evidence was re- 
vealed through observations of X-rays and TeV g amma rays 
from several shell-type SNRs. Kovama et al.l (1 19951) discov- 
ered synchrotron X-rays from the shell of SN 1006, which 
indicates electrons are accelerated up to multi-TeV energies. 
This finding was followed by detections of synchrotron X- 
rays from other SNRs, including RX J1713.7-3946 (e.g. 
iKoyama et all 119971: ISlane et al.l EzOOll) . Further evidence 
for multi-TeV particles (electrons and/or protons) has been 
provided by discovery of TeV gamma rays from so me 
SNRs, such as Cassiopeia A (lAharonian et al.l 1200 ll) or 
RX J 1 7 1 3 . 7-3946 dMuraishi et alJ2000l) . although their spec- 
tral parameters and morphologies were not well determined 
due to the limited sensitivity of TeV observatories. Sub- 
sequently, high quality morphological and spectral studies 
have been p erformed by H.E.S.S. (e.g. lAharonian ^t~aT1l2004l 
l2006ll2007l) . These pioneering measurements by the H.E.S.S. 
telescope, together with the high resolution X-ray data, have 
enabled direct comparison of X-ray and TeV gamma-ray data. 

The shell-type SNR RX J17 13.7-3946 (also known as 
G347.3-0.5), is one of the best-studied SNRs from which 
both non-thermal X-rays and TeV gamma rays are de- 
tected. This SNR was di scovered in soft X-rays during th e 
ROSAT All-Sky Survey dPfeffermann & Aschenbachlll996l) . 
The ASCA satellite, with wider energy coverage than that 
of ROSAT, revealed that the X-ray spectrum is featureless 
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and can be best interpreted as synchrotron emission from 
very high energy electro ns in the TeV regime (Kova ma et alJ 
119971; ISlane et all [l 9991) . The X-ray spectrum was well fit- 
ted with a power-law function of photon index T = 2.2- 
2.4 and interstellar absorption column density Nh = 0.6- 
0.8 x 10 22 cm -2 without any observable evidence for a ther- 
mal emission component. Subsequent observations by Chan- 
dra and XMM-Newton have unveiled structure with a com- 
plex network of bright filaments and knots, in the western 
part of the SNR dUchivama et al.ll2003l; lLazendic et alj|2004t 
ICassam-Chenai et al.ll2004tlHiraga et alj|2005l) . 

TeV gamma-ray emission from RX J17 13.7-3946 was 
first reported by the CANGAROO collaboration in 1998 
(Mura ishi et al.1 uOOOl) . and confirmed by the subsequent 
observations with C ANGAROO-II in 2000 and 2001 
(lEnomoto et al.ll2002l) . Later, the H.E.S.S. collaboration ob- 
tained a resolved im age of the source in TeV gamma rays 
dAharonian et al.l2004l) showing that the gamma-ray emission 
from RX J1713. 7-3946 arises mainly in the shell. These ob- 
servations revealed a striking correlation between the X-ray 
and the gamma-ray images, which indicates a strong con- 
nection between the physical processes resp onsible for X-ray 
and T eV gamma-ray emission components (lAharonian et al.1 
2006). Based on the spectral and morphological information, 
they discussed two possible gamma-ray emission scenarios, 
one where gamma rays are generated by inverse Compton 
scattering of accelerated electrons with diffuse radiation fields 
(the so-called leptonic scenario) and the other where the de- 
cay of secondary 7r°-mesons is responsible for gamma rays 
(hadronic scenario). The later observations with H.E.S.S. re- 
vealed that the flux extends to 30 TeV and, likely, beyond, 
which implies particle acc eleration up to energie s well above 
100 TeV for either model (lAharonian et alj2007l) . 

Most recently, our X-ray observations using Chandra and 
Suzaku have provided important clues for understanding the 
acceleration process in the SNR. From a series of observations 
of the northwest part of the SNR with Chandra in 2000, 2005 
and 2006, we discovered that compact regions of the north- 
west (NW) shell are v ariable in flux on a one-year time scale 
dUchivama et al.l2007l) . The fast variability was interpreted as 
one-year scale acceleration and synchrotron cooling of elec- 
trons with amplified magnetic fields of order of 1 mG. Such 
a large magnetic field in compact regions strongly favors 7r°- 
decay emission as the origin of TeV gamma rays. Also, thanks 
to the wide-band coverage of Suzaku and its low background 
level, we were able to measure a hard X-ray spectrum up to 
40 keV from the southwest portion of RX J17 13.7-3946 with 
a clear ind ication of a high-en ergy cutoff in the synchrotron 
spectrum ( Taka hashi et al.l2008l hereafter Paper I). Combined 
with the upper limit on a shock speed of 4500 km s -1 placed 
by Chandra, the cutoff energy determined by the Suzaku ob- 
servation of the southwest part indicates that particle acceler- 
ation within the SNR shock is so efficient that it approaches 
the theoretical limit corresponding to the s o-called Bohm dif- 
fusion regime (e.g. iMalkov & Drurvll200Tl) . 

In this paper, we present results of mapping observations 
of RX J17 13.7-3946 with Suzaku, which covers about two- 
thirds of the SNR region with 1 1 pointings. The low back- 
ground level of the Hard X-ray Detector (HXD) enables us 
to detect hard X-ray emission up to ~ 40 keV from each of 
the pointings. At the same time, its small field-of-view (FoV) 
of ~ 25' x 25' FWHM gives us information about the spa- 
tial distribution of hard X-ray emission and spectral differ- 
ences from region to region. Thanks to its low instrumental 



background and large effective area, the other detector sys- 
tem aboard Suzaku, the X-ray Imaging Spectrometer (XIS), 
also uncovers new observational facts, such as spectral fea- 
tures below 10 keV and the morphology of relatively dim 
regions left unclear in previous studies by ASCA, Chandra, 
and XMM-Newton. By combining the XIS and HXD spectra 
summed over the data from all the pointings, we show a wide- 
band X-ray spectrum (0.4-40 keV) with quite high statistics, 
with which we investigate not only the existence of a cutoff, 
but also its shape. We then compare the cutoff shape obtained 
with theoretical predictions. 

In Section 2, we describe our Suzaku observations and the 
data reduction procedures. Analysis and results of HXD 
and XIS data are shown in §3.1 and §3.2, respectively. We 
present the wide-band spectrum by connecting the XIS and 
the HXD data in §3.3. A detailed study regarding the cutoff 
structure is also given there. Section 4 is devoted to multi- 
wavelength spectral and morphological studies. The results 
obtained are discussed in the following section, and the re- 
sults are finally summarized. Throughout this paper, we as- 
sume that the d istance to RX J17 13.7 -3946 is close to 1 kpc 
as proposed by Kovama et al.1 (1 19971) based on the Nu value. 
A similar distance has been claimed based on the NANTEN 
CO data dFukui et alj|2003t iMoriguchi et al.1 120051) . The typ- 
ical age of the remnant for such a distance is estimated of 
order of 1000 yr, which can be an indication of association of 
R X J 1 7 1 3 .7—3946 wi th an explosion in A.D. 393 as proposed 
bv lWangetai1dl997l) . 

2. SUZAKU OBSERVATIONS AND DATA REDUCTION 

The Suzaku observatory ( Mitsu da et al.l [2007b is the fifth 
Japanese X-ray astronomy satellite, jointly developed by 
Japan and the US. Its scientific pa y load consists of two co- 
aligned det ector systems, the XI S (Kovama et al. 2007|) and 
the HXD dTakahashi et alj |2007; Koku bun et alj|2007l) . The 
XIS consists of four X-ray CCD camer as which are lo- 
cated in the foci of X-ray telescopes (XRT; ISerlemitsos et al.l 
[2007b . Three of the XIS sensors are front-illuminated (FI; 
0.4-12 keV) CCDs and the other is back-illuminated (BI 0.2- 
12 keV). The non-imaging, collimated hard X-ray instrument, 
the HXD, covers the 10-600 keV bandpass. Two main detec- 
tor units, silicon PIN diodes and GSO scintillators, are buried 
at the bottom of well-type active shield of BGO. The former 
covers the lower energy band of 10-60 keV, while the latter 
detects higher energy photons of 40-600 keV. 

We performed 11 pointing observations of 
RX J1713.7-3946 with Suzaku. The observation log is 
summarized in Table 1, and the pointing position of each 
observation is shown in Figure [TJ The southwest part of the 
SNR, labeled as Pointing 0, was observed in 2005 during 
the Performance Verification phase, while the other ten 
observations were performed in 2006 during the Suzaku AOl 
phase. Since the SNR is located on the Galacic plane, it is of 
importance to check the hard X-ray background associated 
with the Milky Way. We therefore observed two nearby back- 
ground regions containing no bright X-ray point sources. The 
pointing positions of these "OFF" observations are shown 
with red squares in Figure QJ The XIS was operated in the 
normal full-frame clocking mode without spaced-row charge 
injection during all the observations. Since the results from 
Pointing 0, together with those from the OFF observations, 
are already reported in Paper I, we do not give a detailed 
description on the analysis and results for these data. 

We used data products from the pipeline processing version 
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FIG. 1.— Suzaku FoV of each observation of the RX J 1713.7— 1713 re- 
gion o verlaid on the ASCA GIS image (1-5 keV) taken from Uchfv ama et all 
12005). The small squares corresponds to the FoV of the XIS. the outer dia- 
monds drawn for each "OFF" pointing is the rough shape of the FOV (50% 
of the effective area) of the HXD PIN. The numbers indicated in the XIS FoV 
are pointing IDs used throughout this paper. 

1.2. For the XIS analysis, we retrieved "cleaned event files" 
which are screened using standard event selection criteria. For 
the 2006 data, we recalculated the values of pulse invariant 
(PI) and the grade values since incorrect CALDB is applied 
to the pipeline processing of these data as announced by the 
Suzaku instrument teams. We further screened the cleaned 
events with following criteria as recommended by the Suzaku 
instrument teams - (1) cut-off rigidity larger than 6 GV and 
(2) elevation angle from the Earth rim larger than 10°. For the 
HXD data, "uncleaned event files" were screened using stan- 
dard event screening criteria. The exposure times after these 
screenings are shown in Table 1. Due to unstable operation 
of 16 PIN diodes installed in Well-counter units, W00-W03; 
hereafter W0, the bias voltage for these diodes was reduced 
to 400 V from the nominal voltage of 500 V on 26 May 2006. 
On 2006 October 4, the bias voltage of 16 more PIN diodes 
(the PIN diodes in the Well-counter units of W10-W13; Wl) 
was reduced to 400 V for the same reason. The reduction of 
the PIN diode bias voltage leads to a decrease of their effective 
area and also affects their energy response. Since the current 
response matrices do not include these effects, only the PIN 
diodes with a bias voltage of 500 V are utilized in the follow- 
ing analysis. Throughout this paper, the data reduction and 
analysis are performed using HEADAS 6.2 and the spectral 
fitting is done with XSPEC 1 1.3.2. 

3. ANALYSIS AND RESULTS 

3.1. HXD Data Analysis 

3.1.1. Spectral Analysis 

The HXD PIN spectrum from each pointing was con- 
structed and compared with the background model estimated 
for the each observation period. In t he analysis below, the 
non-Xray background (NXB) model dWatanabe et al.ll2007l) 
provided by the HXD team is used for the background genera- 
tion. Since the NXB model does not include the contributions 
from the cosmic X-ray background (CXB), a simulated CXB 



spectrum was added to the NXB model. Specifically, based 
on the r eanalysis of the da ta from HEAO-1 observations in the 
1970's dGruber et al.ll 19991) . the CXB spectrum was modeled 

as 

= 7.9 e ke v- L29 exp (-^-) ph g -i ^V" 1 cm" 2 str^D 
as \ s p J 

where e^y = e/l keV and e p = 41.13. We es- 
timated the CXB signal in each HXD PIN spec- 
trum using the response matrix for spatially uniform 
emission, ae_hxd_pinflat_20060809 . rsp, 

ae_hxd_pinflatl23_200 6080 9 . rsp, and 
ae_hxd_pinflat23_20060809 . rsp. The contri- 
bution from the CXB flux to the detected count rate is 
estimated to be — 5% of the NXB. 

Figure [2] shows the HXD PIN spectrum obtained from 
Pointing 8, where a clear detection of hard X-rays can be 
seen. Likewise, the HXD PIN detected signals from all other 
pointings with a count rate of 20-50% of the NXB. We fitted 
the background-subtracted spectra with a simple power law 
model: dN/de oc e~ r . Since the extended nature of the source 
does not cause any spectral steepening or flattening (Paper 
I), we used the point-source response matrix at the XIS- 
nominal position, ae_hxd_pinxinom_20060814 .rsp, 
ae_hxd_pinxinoml23_2 00 60814 . rsp, and 
ae_hxd_pinxinom23_2 00 60814 . rsp. Table 2 
gives the best-fit parameters with the statistical errors at 90% 
confidence level. The obtained photon indices are generally 
larger than those obtained from the corresponding XIS 
spectra (see below). This difference indicates that a spectral 
cutoff is not unique to the S W region (Paper I), but a common 
feature throughout the remnant. 

In order to confirm the results obtained above, we evalu- 
ated the systematic errors due to uncert a inties in the NXB 
modeling. As described in iMizuno et al.l (12006). the current 
reproducibility of the NXB model is ~ 5%. Therefore, we ex- 
amined how much the values of the photon index change by 
increasing or decreasing the background model by 5%. The 
systematic errors were found to be smaller than the statistical 
errors indicated in Table 2. 

The systematic errors due to the misestimation of NXB 
were examined in another way. Considering the physical 
size of this target, the emission from the remnant should be 
constant during the observations. Therefore, background- 
subtracted lightcurves should be constant during a observa- 
tion. Although the lightcurves shown are almost constant 
within statistical errors, the background-subtracted count rate 
becomes higher when the total count rate increases for Point- 
ings and 8. The light curves for Pointings and 8 are 
shown in Paper I and Figure [3] respectively. Since this be- 
havior is thought to be caused by misestimation of the NXB, 
we examined how much the photon indices change with and 
without those time regions. When we discard the time region 
which corresponds to the last bin of Figure[3] the photon index 
changes by AT ~ 0.2 from the values in Table 2. 

Diffuse emission from the Galactic plane also can affect 
the spectra. However, no emission above the 5% level of the 
NXB was detected from the OFF pointings (Paper I). More- 
over, when the excess counts marginally detected from the 
OFF pointings are added to the background spectrum for each 
pointing, the fitting result agrees within the statistical errors. 

3.1.2. Spatial Distribution of Hard X-ray Emission 
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TABLE 1 

Summary of the Suzaku observations of RX J 1 7 1 3 .7— 3946 



Pointing ID 


Obs. ID 


Coord. (J2000) 


Exposure 


Date 






("J2000. <5j2000) 


XIS/HXD 










[ks] 








100026010 


(17 h 12 m 17 s 


.0. 


-39 d 56 m ll s ) 


55/48 


26/9/2005 


1 


501063010 


(17 h ll m 51 


.5, 


-39 d 31 m 13 s ) 


17/17 


11/9/2006 


2 


501064010 


(17 h 12 m 38 s 


■0, 


-39 d 40 m 14 s ) 


18/22 


1 1/9/2006 


3 


501065010 


(17 h 12 m 38 s 


•2, 


-39 d 22 m 15 s ) 


19/18 


1 1/9/2006 


4 


501066010 


(17 h ll m 04 ! 


.5, 


-39 d 40 m 10 s ) 


19/21 


12/9/2006 


5 


501067010 


(17 h ll m 05- 


■1, 


-39 d 22 m 10 s ) 


16/19 


12/9/2006 


6 


501068010 


(17 h 14 m ll s 


.6, 


-39 d 40 m 14 s ) 


20/19 


13/9/2006 


7 


501069010 


(17 h 14 m ll s 


•4, 


-39 d 22 m 15 s ) 


12/11 


19/9/2006 


8 


501070010 


(17 h 14 m ll s 


.8, 


-39 d 58 m 14 s ) 


19/20 


19/9/2006 


9 


501071010 


(17 h 12 m 17 s 


.6, 


-39 d 18 m 50 s ) 


16/15 


20/9/2006 


10 


501072010 


(17 h 15 m 44- 


.5. 


-39 d 40 m 10 s ) 


15/15 


5/10/2006 


OFF1 


100026020 


(17 h 09 m 31 


.9, 


-38 d 49 m 24 s ) 


28/24 


25/9/2005 


OFF2 


100026030 


(17 h 09 m 05 


• 1, 


-41 d 02 m 07 s ) 


30/28 


28/9/2005 



TABLE 2 

Power-law fitting to the HXD PIN spectra" 



Pointing ID 


r 


F[(h40 keV 
(mCrab) 




Used Units 





3.2 ±0.2 


2.5±0.1 


1.15 (36) 


W0-3 


1 


3.3 ±0.2 


3.2±0.2 


0.91 (18) 


Wl-3 


2 


3.0 ±0.3 


3.2 ±0.2 


0.84 (21) 


Wl-3 


3 


3.4±0.5 


2.0 ±0.2 


0.58 (17) 


Wl-3 


4 


2.9 ±0.3 


3.5 ±0.2 


0.77 (23) 


Wl-3 


5 


3.2±0.5 


1.7 ±0.2 


0.76 (17) 


Wl-3 


6 


2.9 ±0.3 


2.9 ±0.2 


0.72 (20) 


Wl-3 


7 


3 q+1.4 
"-LI 


1 2+ 03 
'■^-0.2 


0.42 (9) 


Wl-3 


8 


2 6 +a5 


2.2 ±0.2 


0.87 (19) 


Wl-3 


9 


3.0+ 1 - 2 


1.1 ±0.2 


0.86 (12) 


Wl-3 


10 


A 4+1 : <> 


1.0±0.2 


0.20 (8) 


W2, 3 



' Errors represent 90% confidence. 




20 30 
Energy [keV] 

FIG. 2. — The HXD PIN spectrum obtained from the observation of Point- 
ing 8. The data points in black show the raw spectrum, red points represent 
the background model (NXB+CXB), and green points are the background- 
subtracted spectrum. 

Using the spectral parameters obtained for the 11, pointings 
we attempt to reconstruct the spatial distribution of the hard 
X-ray emission. Since we expect spatial variation of both 
brightness and spectral index, we use a Monte Carlo simu- 
lator, simHXD (Ter ada et al.l 12005). for the estimation. We 
input the spatial distribution of brightness and spectral index 
into simHXD, and compare the simulated spectra with ob- 
served ones shown in 33.1.11 



r 0.4 




10000 20000 30000 40000 

Time [s] 

FIG. 3. — The HXD PIN light curve during the observation of Pointing 8. 
Events with energy from 12 keV to 40 keV are selected. The black, red, and 
green points show the raw, background (NXB), and background-subtracted 
data, respectively. 

Firstly we simulated by taking a simple model in which the 
brightness distribution is same as the ASCA image shown in 
Figure|4]and the spectral index is constant at T = 3.0 through- 
out the SNR (simulation!). By using the ASCA image as an 
input to simHXD, we simulated hard X-ray spectra in 10-40 
keV, which are expected to be observed by the PIN in each 
pointing. Figure [5] compares the flux from the observations 
and the simulations. In this figure, each value is normalized 
to that of Pointing 0. Here, the systematic error for each ob- 
servational data point is ~ 20% if we use the 5% of the back- 
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FIG. 4. — The emission distribution (left) and the photon index distribution 
(right) adopted in the simulations. The emission distribution is the ASCA 
image in the energy range of 1-5 keV. 

ground as systematic errors. The observations and the simula- 
tions appear consistent with each other. Therefore, we expect 
the brightness distribution above 10 keV is not drastically dif- 
ferent from the distribution below 10 keV. 

The largest discrepancy between the observation and the 
simulation is found for Pointing 4. In this pointing, 
there is a known point source at the corner of the HXD- 
PIN field of view. The source is listed in t he Second 
IBIS/ ISGRI Soft Gamma-Ray Survey C atalog (iBirdetalJ 
[2006b as IGR J17088-4008. According to lBird et all 12006), 
the average flux of this source is 1.1 ±0.2 mCrab in 20- 
40 keV and 2.2 ± 0.3 mCrab in 40-100 keV. It is noted that 
the source is not bright enough for IBIS/ISGRI to determine 
spectral parameters. The transmission of fine collimators of 
the HXD for the source is estimated to be ~ 0.05. The esti- 
mated count rate of this source corresponds to ~ 2% of the 
detected signals from Pointing 4. However, the variability of 
this source could increase the count rate to a non-negligible 
level. The observation and the simulation become consistent 
if the point source was 10 times brighter than its average value 
during the Suzaku observation. Also, the angular response of 
the HXD near the edge of the FoV can have some uncertain- 
ties since calibrations for a source with such a large offset 
angle is difficult. The uncertainties can affect our estimate of 
the contamination from IGR J 1708 8-4008. At this moment, 
we cannot conclude that the large difference between the ob- 
servation and the simulation is due to the spatial distribution 
of the SNR emission or the point source in the FoV. 

Next we tried another simulation taking into account the 
spatial distribution of spectral indices (simulation2). We 
adopt a "toy model" shown in Figure [4] since the results of 
the spectral analysis presented in Table 2 suggest that the hard 
X-ray spectrum may be flatter in the inner region of the SNR 
than near the rim. In the model, the photon index of the in- 
ner region is set to T = 2.6 and that of the rim region is set to 
r = 3.5. The ASCA image was used to provide the brightness 
distribution at 10 keV. 

Comparisons of the data and the simulation results are 
shown in Figure [5] and Figure [6] for the flux and photon index, 
respectively. The detected flux obtained from simulation! is 
somewhat similar to that from simulationl , and the simula- 
tion data generally follow the observational data. As for the 
distribution of photon indices, the toy model gives a similar 
distribution to the observational results for the western por- 
tion of the SNR. Obtaining a better fit for the eastern region 
(Pointing 7, 8, and 10) may require more complex assump- 
tions than the toy model for simulationl. 

3.2. XIS Data Analysis 
3.2.1. Image Analysis 
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FIG. 5. — Comparison of detected flux obtained from the observations, 
simulationl, and simulation!. Each data point is normalized to the value of 
Pointing 0. The error bars indicate la statistical errors. 
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FIG. 6. — Comparison of photon indices obtained from the observations 
and simulation!. The error bars correspond to lcr statistical errors. 

Figure [7] shows mosaic images of RX J1713. 7-3946, con- 
structed using the data from XIS0, 2, and 3 (FI-CCDs). The 
upper panel shows the soft band image in 1-5 keV and the 
lower panel shows the hard band image in 5-10 keV. Both 
images are smoothed with a Gaussian of a = 0'.3. Instru- 
mental background signals are subtracted from both images. 
The signal to background ratio in the hard band is smaller 
than that in the soft band by one order of magnitude. Thus, 
the background must be carefully subtracted. We utilized the 
Night Earth Background Database consisting of event data 
obtained when the satellite is looking at the night earth and the 
non X-ray background becomes dominant. After subtracting 
the background, the vignetting effects of the XRTs were cor- 
rected by means of the simulation program called xissim 
(Ishisa ki et alJ 120 07). In the program, an image from a flat 
field can be simulated by a Monte Carlo method. 

As is seen in Figure [7] the Suzaku XIS has covered most 
of the remnant. Thanks to little stray-light contamination of 
the XRTs and low background level of the XIS, high qual- 
ity images are obtained even in the high energy band above 
5 keV. The double shell structure revealed by XMM-Newton is 
clearly seen in the XIS images. In addition to the bright struc- 
tures of the western part, the XIS revealed detailed morphol- 
ogy of the dim parts of the remnant. The dim structures are 
highlighted in Figure [8] which is the same as Figure [7] (a) but 
displayed with a different color scale. In thi s figure, the con- 
tours of the H.E.S.S. gamma-ray image (from lAharonian et al.l 
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120071) are overlaid for comparison. As is clearly seen, not only 
the bright rims but also the eastern portion shows a striking 
similarity between the two energy regimes. This correlation 
is discussed more quantitatively in 33.41 

Two point sources seen in the soft band image are 
listed in the ROSAT bright source catalogue. One is 
1WGA J171 44-3945, which was associated w ith a Wolf- 
Rayet star by Pfeffermann & Aschenbach ( 1996). The other 
is 1WGA J17 13.4-3949, which is located between the two 
FoVs of Pointing 2 and 6. This source has be en suggested to 
be the neutron star associated with the SNR by Lazendic et al. 
J2Q031 

Comparing the soft-band and hard-band images provides 
information about the spatial variations of the spectral proper- 
ties. At first sight, the bright structures are very similar in the 
soft band and hard band images. In order to compare the im- 
ages in more detail, a radial profile around the center of SNR 
(«j20oo = 17 h 13 m 33 s .6, feooo = -39 d 45 m 36 s ) is presented in 
Figure|9] where circular regions of 2' . 1 radius centered on the 
two point sources are excluded. The hard band to soft band 
ratio is significantly different between the bright outer region 
part and the interior. This difference suggests a corresponding 
difference in spectral properties. 

3.2.2. Spectral Analysis 

As already reported in Paper I, the XIS data reveal a cut- 
off below 10 keV in the spectrum of the SW rim. The XIS 
spectra from the various regions should reveal whether the 
cutoff varies across the remnant. We extracted spectra from 
the regions shown in Figure [7] with the green squares. The 
region labeled as "1' (Region 1) is the same one that we 
used for the spectral analysis in Paper I. The background 
spectrum needs to be accumulated from nearby regions since 
RX J1713. 7-3946 is located almost on the Galactic plane. 
The OFF observation data cannot be used for background es- 
timation other than for Region 1 . This is becaus e the contami- 
nation on the optical blocking filters of the XIS ( Koyama et al. 
2007) significantly changed the detector response to low- 
energy X-rays between 2005 and 2006. For Region 1-10 we 
therefore accumulated a background spectrum from the re- 
gions indicated with the magenta polygons in Figure [7] In 
the spectral fitting discussed below, the standard RMF files 
version 2006-02-13 were used, whereas ARF files were pro- 
duced using xissimarf gen. Each spectrum is binned so 
that each bin contains at least 300 counts. After the binning, 
we ignore those bins whose energy is smaller than 0.4 keV or 
larger than 12 keV. We also excluded bins between 1.7 and 
1.9 keV because there exist large systematic uncertainties in 
the response matrices. Finally, we co-added the spectra, RMF 
files, and ARF files from the FI chips to produce a single data 
set for the XIS. 

First, we fitted all the spectra with a simple power law ab- 
sorbed by the interstellar medium. The results are summa- 
rized in Table 3. This model yields acceptable fits for most 
of the spectra. However, the fits to the spectra of Region 1, 
2, and 6 are not acceptable even at the 99% confidence level. 
We show the spectrum of Region 2 with the best-fit power- 
law and the residuals in Figure [10] Although not so drastic 
as Region 1 which is shown in Paper I, a correlated pattern to 
the residuals can be seen for the Region 2. A similar pattern 
is also seen in the data of Region 6. 

We the fitted all the spectra with a cutoff power law. This 
function gives acceptable fits for all regions, and consistently 
gives better values of xl- In Figure [10] we plot the residuals 
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FIG. 7.— The XIS (XIS 0+2+3) mosaic images of RX J1713.7-3946 in 
the energy band 1-5 keV (a) and 5-10 keV (b). North is up and east is to 
the left. The color scale indicates count rate on a linear scale. The green 
regions in the top panel are those used for spectral analysis. The dashed 
magenta polygons correspond to the regions used to extract the background 
spectrum. The positions of the two point sources, 1WGA J1714.4— 3945 and 
IWGA J1713.4— 3949 are shown with white circles. The dashed green cir- 
cles represent radii of 0.3° and 0.6° around the center of the SNR (aj2000 = 
17 h 13 m 33.6 s , <5j2ooo = -39 d 45 m 36 s ), which correspond to the vertical dashed 
lines in Figure|9] 

for the spectrum of Region 2. These results, together with the 
steeper spectra detected with HXD, indicate the existence of 
a cutoff somewhere between the bandpasses of the XIS and 
HXD. 

Although statistically rejected for some regions, the re- 
sults obtained with the power-law fits are consistent with 
previous studies wi t h ASCA, Chandra, and XMM-Newton 
dKovama et al.lll997tfsTane et alJl999tlUchivama et al.ll2003b 
ICassam-Chenai et al.ll2004t iHiraga et alj |2005l) . The results 
with Suzaku suggest that the photon index is larger and the 
absorption is smaller for the inner regions than for the outer 
regions. This tendency agrees with the radial profile shown in 
Figure [9] The difference of /Y H between the western bright 
spots and the inner region is A/Vh — 0.3-0.4 x 10 22 cm" 2 , 
which is consistent with the results from XMM- Newton. Ac- 
cording to the discussion by Hirag aetail d2005l) . there could 
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FIG. 8. — Comparison of the Suzaku XIS image and the g amma-ray im - 
age by the H.E.S.S. telescope (contours) taken from Aharoni an et al.l j2007l) 
shown with color scale and green contours, respectively. The XIS image is 
same as Figure UJ (a) but the scale is changed to stress the similarity of the 
two images. 




0.4 

Radius [deg] 

FIG. 9. — Radial profiles around the center of the SNR (aj2000 = 
17 h 13 ra 33.6\ feooo = -39 d 45 m 36 s ) in the two energy bands. Each profile 
is normalized to its peak value for direct comparison. The radii correspond- 
ing the vertical dashed lines are shown in Figure □(b) with the green dashed 
circles. 

be a correlation between the difference of Nh and the pres- 
ence of the molecular clouds in the weste rn part of th e SNR 
detected with the N ANTEN telescope dFukui et all 120031: 
iMoriguchi et"aT1l2005l) . 

3.3. Broad-Band X-Ray Spectral Analysis 

In this section, we connect the XIS and HXD spectra, which 
is crucial for a quantitative study of the cutoff structure and 
also for multi-wavelength study described in §13.41 In or- 
der to study the general characteristics of the emission from 
RX J1713.7-3946, the XIS and HXD spectra obtained in 
CLU'SLDwere co-added. Then, the summed XIS and HXD 
spectra were scaled to account for the flux from the whole 
remnant, which makes it easier not only to connect the XIS 
and HXD spectra but also to compare the combined spectrum 
directly to those of other wavelengths in § 13.41 In scaling the 
spectra to the whole remnant, we assumed the surface bright- 
ness of the ASCA GIS image (1-5 keV) shown in FigureQ] the 
angular response of the XIS/XRT system (ARFs), and that 
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FIG. 10.— The XIS (XIS 0+2+3) spectrum from Region 2. The lower 
panels show the residuals when the spectrum is fitted with a power law and a 
power law with an exponential cutoff. 



of the HXD PIN presented in Figure 6 of Paper I. The rela- 
tive normalization factor between XIS and HXD derived from 
Crab observations (Ishida et al. 2006) is included in the scal- 
ing. Therefore, the scaled XIS and HXD can be connected 
to each other without additional scaling if no systematic er- 
rors are considered. However, we estimate that the normal- 
ization of the XIS and HXD spectra obtained by the proce- 
dures above should contain systematic errors of 10-20%. In 
order to account for the systematic errors, we include a con- 
stant normalization factor and deal with this factor as a free 
parameter when fitting the XIS+HXD spectrum below. In this 
procedure, we discarded the HXD data from Pointing 4 since 
they seem to be contaminated by a nearby hard X-ray source 
as described in §13.1.21 

Before the XIS and HXD spectra were jointly analyzed, 
each co-added spectrum was independently fitted using a 
power-law function. Table 4 summarizes the fit results. The 
HXD spectrum gives an acceptable fit. In contrast, the XIS 
does not, with a large chi-squared value of Xu 2 = 1 -57 for 71 1 
degrees of freedom and residuals that begin to become large at 
~ 6 keV. This fact is clearly seen in FigureQTl where the XIS 
and HXD spectra are plotted together with a power-law func- 
tion which represents the XIS spectrum. The spectral steep- 
ening begins in the XIS band and continues smoothly into the 
HXD band. This plot strongly suggests that spectral steep- 
ening occurs around 10 keV. We have already reported such 
a spectral feature in Paper I for Pointing data. The same 
kind of feature has been revealed in this spectrum averaged 
over the SNR, which suggests that the spectral steepening is 
common in the entire region. 

We then quantitatively evaluated the spectral steepening 
by fitting a model with a cutoff structure to the combined 
XIS/HXD spectrum. Below we use our numerically calcu- 
lated synchrotron spectrum which was embedded into the 
XSPEC package. As an electron distribution, we adopt a gen- 
eralized form, namely, 



dN e 
~d~E 



oc E s exp 



-) 



(2) 



instead of taking a simple exponential cutoff often used in the 
literature. Here, s represents index of electron spectrum and 
/3(> 0) determines rapidity of high-energy cutoff. The photon 



8 



Tanaka et al. 



TABLE 3 

Model fitting to the XIS spectra" 



Region 



Area 



Power law 



Cutoff power law 









r 


^1-10 kev" 






r 




Fl-10 kev" 






(arcmin 2 ) 


(10 22 crrT 2 ) 




(10~" erg s" 1 cm" 2 ) 




(10 22 cm 4 ) 




(keV) 


(10 -11 erg s _1 cm -2 ) 




1 


108 


0.87 ±0.01 


2.39 ±0.01 


6.06 ±0.03 


1.38(718) 


0.77 ±0.01 


1.96 ±0.05 


9±1 


5.7±0.1 


1.04(717) 


2 


77.2 


0.83 ±0.02 


2.38 ±0.02 


5.22 ±0.05 


1.23 (242) 


0.73 ±0.03 


2.00 ±0.09 


10+j 


4.9 ±0.2 


1.01 (241) 


3 


86.6 


0.55 ±0.03 


2.62 ±0.04 


1.41 ±0.03 


1.14(94) 


0.47 ±0.04 


2.2 ±0.2 


8+ T 


1.3 ±0.1 


1.02 (93) 


4 


42.2 


0.54 ±0.04 


2.72 ±0.07 


0.66 ±0.02 


1.38 (44) 


0.44 ±0.07 


2.2 ±0.3 




60 m06 


1.23 (43) 


5 


116 


0.78 ±0.03 


2.39 ±0.04 


1.88 ±0.04 


1.18(120) 


0.70 ±0.05 


2.1 ±0.2 


id 3 


1.8±0.1 


1.12(119) 


6 


51.4 


0.88+0;°! 


2.43 ±0.06 


1.27 ±0.04 


1.63 (48) 


0.68 ±0.09 


1.7±0.3 




1 l 40 - 2 
la -o.i 


1.24(47) 


7 


59.6 


0.71 ±0.04 


2.38 ±0.05 


1.18±0.03 


1.15 (56) 


0.62 ±0.07 


2.0 ±0.2 


10+' 4 


1.1 ±0.1 


1.03 (55) 


8 


72.4 


0.78 ±0.03 


2.55 ±0.04 


1.91 ±0.04 


1.10(110) 


0.68 ±0.05 


2.1 ±0.2 


4 


1.8 ±0.1 


0.92 (109) 


9 


65.9 


0.81 ±0.03 


? 4Q+0.04 


2.03 ±0.04 


1.11 (109) 


74+ 04 


2.2 ±0.2 


j 4+1 6 


1 Q+0.2 

1 - y -o.i 


1.03 (108) 


10 


54.0 


0.60 ±0.04 


2.23 ±0.05 


1.49 ±0.03 


0.88 (75) 


0.52 ±0.06 


1.9±0.2 


12« 9 


1.4±0.1 


0.79 (74) 



1 Errors represent 90% confidence. b Corrected for absorption. 

spectrum can be calculated as 

as J \e c J dE 
Here the function F(x) is denned as 

/•OO 

F(x)=x K 5/3 (Qd£, 



(3) 



(4) 



where is the modified Bessel function of 5/3 order. When 
pitch angles are isotropic, the characteristic photon energy e c 
is given as 



e c = 0.543 



B 



100 J V 10 TeV 



keV. 



(5) 



The model consists of four parameters: s, H = EqB 1 ! 1 , (3, and 
the flux normalization. Photon spectra were calculated and 
tabulated for reasonable ranges of the four parameters, to be 
used as an XSPEC table model. The parameters characteriz- 
ing the electron energy distribution can be obtained directly 
from the fit to the X-ray data. 

For RX J17 13.7-3946, the most probable value for s is 3.0 
rather than 2.0 due to significant synchrotron cooling of elec- 
trons during the lifetime of the SNR to (~ 1000 yr). The elec- 
tron spectrum becomes steeper by a factor of 1 (5 — > 
when the injection of electrons is constant and the lifetime 
to is smaller than the time scale of synchrotron cooling f sync . 



Since f sync can be written as 



( B 



fsync " 28 V100//G 



-3/2 



Gkev) 



-1/2 



yr, 



the condition above is satisfied in the energy range, 



£ > El, = 2.3 



B 



( to 



100 nG) V 103 y r 



eV. 



(6) 



(7) 



To explain the variability of X-ray emission reported on 
a year timescale from compact r egions of the shell of 
RX J17 1 3.7-3946. lUchivama et al.l (120071) proposed that the 
magnetic field in these compact regions is amplified to 1 mG. 
The average large-scale magnetic field in the remnant should 
be significantly lower, but even for B ~ 100 fiG the above 
condition in the X-ray domain is safely satisfied. Therefore, 
for a strong shock with a compression ratio of 4.0, the index 
of X-ray emitting electrons should be close to s = 3.0. 

We fitted the Suzaku spectrum with the synchrotron spec- 
trum described above. In the fitting procedure, we fixed the 



electron index s. Table 5 summarizes the result in the case of 
s = 3.0, which corresponds to the most probable case follow- 
ing the discussion above. We also present results when s = 2.0 
for comparison. This case can be realized if the magnetic field 
and/or the age of the remnant are much smaller than we ex- 
pect, as is seen in equation (0. The spectral fitting with the 
electron index of s = 3.0 yields rather rapid steepening with 
j3 = 3.4^q g compared to the conventionally used exponential 
cutoff. However, one should be careful with the physical in- 
terpretation of the fit based on equation 0. Indeed, while the 
electron spectrum is derived from the broad-band Suzaku X- 
ray data assuming a specific spectral form given by equation 
10, the curve 1 in Figure [T2] has a more general meaning. It 
does not depend on the assumed analytical presentation and, 
in fact, relates uniquely the energy spectrum of electrons to 
the measured X-ray spectrum. Indeed, the curve 1 can be 
presented in different mathematical forms. In particular, the 
electron spectrum shown by curve 1 in Figure[l2]is quite close 
to the theoretical prediction for the spectrum of shock accel- 
erated electrons in a young SNR (Zirakashvili & Aharonianl 
2007), 



dN e 
~dE 



oc E 



1+0.661 !) 



5/2' 



- 9/5 








exp 


[-©1 



(8) 



This spectrum is derived under the assumption that elec- 
trons are accelerated by a strong shock in the Bohm diffusion 
regime and that the energy losses of electrons are dominated 
by synchrotron cooling. It is seen that the energy spectrum 
below the cutoff is described by a function which deviates 
from a pure power-law form, by an exponential term with 
(3 = 2. The latter has a simple physical interpretation and is 
a result of combination of two effects — acceleration in the 
Bohm diffusion regime, and energy losses in the synchrotron 
regime. The curve 2 shown in Figure 12 is calculated for the 
parameter n = £ B 1/2 = 201 TeV ^G 1 / 2 . 

The energy spectrum of synchrotron radiation correspond- 
ing to the electron spectrum give n by equation ([8]) has a form 
( Zira kashvili & Aharonianl2007l) 



dN 

de 



oc e 



1 + 0.46 



so 



0.6 



-1 2.29 




exp 









(9) 



This spectrum fits well the broad-band Suzaku X-ray data for 
a single parameter £0 = 0.67 ± 0.02 keV. Note that the expo- 
nential term in the synchrotron spectrum is a weak function 
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Fig. 11.— The Suzaku QQS+HXD) spectrum of RX J1713.7-3946. The 
model plotted with the data is a power law obtained by fitting the XIS data. 
The parameters are shown in Table 4. The lower panel shows the ratio of the 
data to the model, and clearly reveals a cutoff. 




E(B/100 mG) 



FIG . 12. — Comparison of electron spectra in the high-energy cutoff region. 
Curve 1 is the electron spectrum obtained by fitting the Suzaku spectrum with 
a synchrotron spectrum assuming the electron distribution given by equation 
[2] whe n s is fixed to 3.0. Curve 2 is the model by Zirakashvili & Aharonian 
12007). Both the spectra are drawn by assuming the best-fit parameters of 
the Suzaku spectrum, which are shown in Table 5. Curve 3 is the exponential 
cutoff spectrum often assumed in the literature. The cutoff energy is set to the 
same value as that of curve 1 . 

of energy (oc exp[-(s/eo) 1 ^ 2 ])', therefore the characteristic en- 
ergy £0 only formally can be considered as the cutoff energy. 
In fact, for the shape given by equation J9), the break in the 
spectrum starts at much higher energies (e ~ lOeo). 

3.4. Multi-Wavelength Study 
3.4. 1. Spectral Energy Distribution 

The high quality data of Suzaku obtained over two decades 
in energy and combined with the TeV gamma ray data, allow 
definite conclusions concerning the origin of multi-TeV par- 
ent particles based on the comparison of model predictions 
with observations. In particular it is important that the Suzaku 
data provide unambiguous information about the shape of the 
energy spectrum of electrons in the cutoff region. 

Figure Qj] shows the Spectral Energy Distribution (SED), 
E 2 dN/dE, of RX J17 13.7-3946 from radio to TeV ener- 
gies. The X-ray data correspond to the entire remnant; they 
are reconstructed assuming the best-fit model given by equa- 
tion (O corrected for the interstellar absorption with Nh = 
0.70 x 10 22 cm" 2 . The points in the TeV gamma-ray range 
corresponding to the fluxes of the whole rem nant, are from the 
latest report of the H.E.S.S. collaboration (Aharonian et al 
2007). The EGRET upper limit is taken from Aharonian et al 
(2006); it is obtained through modeling and subtracting 



the re solved EGRET source 3EG 1714-3857 dHartman et alj 
1999). The two data points in the ra dio band are shown a s 
measured with the ATCA telescope dLazendic et al.1 120041) . 
But it should be noted that these fluxes are detected only 
from the NW rim region, therefore they should be treated 
as lower limits when compared to model predictions for the 
whole remnant. The recent estimates based on observations 
with ATCA and the 30-m radiotelescope IRA show that the 
flux of the entire remnant should be a factor of few larger than 
the flux of the NW rim (F. Acero et al. 2008, in preparation; 
G. Dubner 2008, private communications). 

In Figure[T3]we show theoretical fluxes of broad-band elec- 
tromagnetic radiation produced by both accelerated electrons 
and protons. The results depend strongly on the strength of 
the average magnetic field B. Although the magnetic field in 
comp act filaments can be as large as ~ 1 mG (Uchiya ma et al.1 
120071) . the field in less bright and more diffuse regions where 
the bulk of the synchrotron X-ray emission is produced, 
should be significantly weaker. Yet, a moderately large field 
exceeding 100 /iG is required to explain the X-ray flux ratio 
between the diffuse and compact zones. The integrated X- 
ray flux from the diffuse zone (zone 1) is larger by a factor 
of ~ 4 than that from the compact zone (zone 2): F\ ~ AFi 
dUchivama et al. 2003). Using the ratio of the emission vol- 
umes estimated from the X-ray image (Vi ~ IOOOV2), the 
number density of X-ray emitting electrons, n oc F/(VB 2 ), 
in the zone 2 relative to the zone 1 is estimated as «2 ~ 
250(Bi/B2) 2 «i- For B > 1 mG, and assuming m > «i, one 
finds Bi > 100 /iG. 

In Figure [13] the synchrotron, inverse Compton (IC), and 
7r°-decay fluxes are calculated assuming strong uniform mag- 
netic field of strength B = 200 /iG. The calculations are per- 
formed with constant injection of electrons and protons over 
the last 1000 yr. The injection spectrum for both electrons 
and protons is assumed to be a power-law with an index 
s = 2.0. For calculation s of IC fluxes we use d the diffuse pho- 
ton fields proposed by iPorter et all (120061) . including 2.7 K 
CMB, as well as interstellar radiation consisting of two, opti- 
cal (starlight) and infrared (dust) emission components. 

The best fits for the synchrotron radiation and 7r°-decay 
gamma rays are achieved for cutoffs in the energy distribu- 
tions of electrons and protons at 28.4 TeV and 130 TeV, re- 
spectively. The spectra of 7r°-decay gamma rays are cal- 
culated using the analytical presentations by iKelner et al.1 
d2007l) . Since the highest energy tail of the observed gamma- 
ray spectrum is best fi tted with exp[-(e/eo) t3n >] with /3 7 w 0.5 
(Aharonian et al. 2007), the proton spectrum in the cutoff re- 
gion ^shojfldjiave anexponential behavior, (3 P ~ 2 x (3~ = 1 
(see lKelneret~aTj|2007l) . 

The electrons suffer synchrotron losses, thus the electron 
spectrum above a certain energy becomes steeper (E~ s — » 
fjjg position of the break in the electron spectrum 
appeai-s at 

Eb = 1-25 ( -r^-r; ) (t:^—) TeV. (10) 



100 nGj V 103 y r / 

For a magnetic field B = 200/iG and age of the source T < 
1000 yr, the spectral break in the synchrotron spectrum cor- 
responding to the transition of the electrons spectrum from 
uncooled to cooled regime, appears around 1 eV (see Fig- 
ure [T3l . Thus, for young SNRs the detection of the syn- 
chrotron break at optical/infrared wavelengths would be an 
additional argument in favor of strong magnetic field. The 
cooling break in the electron spectrum is reflected also in the 
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TABLE 4 

Power-law fitting to the XIS/HXD spectrum of the whole remnant 11 



Detector 




r 


Flux b 


x 2 (y) 




(10 22 cm- 2 ) 




(lO^ 10 ergs" 1 cm" 2 ) 




XIS (0.4-12 keV) 


0.79 ±0.01 


2.39±0.01 


7.65 ±0.03 


1.57 (711) 


HXD (12^40 keV) 




3.2±0.1 


1.17±0.03 


1.16(67) 



a Errors represent 90% confidence. 11 Corrected for absorption. The calculated energy 
band is 1-10 keV and 10-40 keV for XIS and HXD, respectively. 



TABLE 5 

Fitting to the Suzaku (XIS+HXD) spectrum 1 



Assumed Function N\\ Constant Factor Parameter(s) ^l-l0keV a Xv 2 (*0 
(10 22 cm -2 ) (lO" 10 erg s~ ] cm" 2 ) 

Electron Distribution 0.71±0.01 1.03 ±0.06 s = 3.0 (fixed), (3 = 3.4+°], n = 402 ± 6 TeV tiG 1 / 2 7.2±0.1 1.11 (778) 

of Equation l|2) c 0.68 ±0.01 1-03_qq5 s = 2.0 (fixed), g = 1.5 ±0.2, n = 207+ 2 ^ TeV /jG 1 / 2 7.1 =1=0.1 1.07 (778) 

Z&A(2007) d 0.70±0.01 1.08 ±6.04 e = 0.67 ± 0.02 keV " 7.2±0.1 1.11 (779) 

a Errors represent 90% confidence. ^ Corrected for absorption. c Details of the calculation of synchrotron spectrum is described in the text. d A model by 
Zirakashvili & Aharonian 1 2007). the synchrotron spectrum of which is shown as equation (9) in this paper 



gamma-ray band. Namely, for B = 200/iG the corresponding 
signature in the IC gamma-ray spectrum appears around 10 
GeV. Unfortunately (see Figure [T3T l. for such a strong mag- 
netic field, the IC component is suppressed, and falls well 
below the sensitivity of gamma-ray detectors in this energy 
band, including GLAST. At the presence of such a strong 
magnetic field, the only viable mechanism which can pro- 
duce TeV gamma-rays at the flux level detected by H.E.S.S. 
is related to interactions of ultrarelativistic protons with am- 
bient gas through production and decay of 7r°-mesons. To 
explain the detected flux of gamma-rays protons should be 
accelerated to energies well beyond 100 TeV, and the parame- 
ter A = (W p /10 50 erg)(n/l cm" 3 )(d/lkpc)" 2 should be between 
1.5-3, depending on the spectrum of protons. 

In order to increase the IC flux to the level of the observed 
TeV gamma-ray flux, the magnetic field should be reduced 
down to 10 fiG and 15 /iG. For the given magnetic field, 
the high quality X-ray data of Suzaku obtained over two en- 
ergy decades allow derivation of the electron spectrum with 
high accuracy within the interval covering one energy decade: 
from - 50(5/10 /xG)" 1 / 2 TeV to - 500(B/10 /xG)" 1 / 2 TeV. 
This allows us to calculate the spectrum and absolute flux of 
IC gamma rays above a few TeV without any model assump- 
tions, as long as the main target for the IC gamma-ray produc- 
tion remains the 2.7 K CMB. The contribution of the diffuse 
optical/infrared radiation fields generally is less, however the 
optical photons may provide enhanced TeV emission at low, 
sub-TeV energies. For calculations of the IC spect r um w e 
used the interstellar radiation model of iPorter et al.1 (120061) . 
who proposed significantly larger flux of optical and infrared 
components compared to the generally accepted flux. How- 
ever, the results presented in Figure [14] show that even this 
high diffuse optical and infrared radiation fails to account for 
the observed gamma-ray flux below a few TeV. In order to fill 
this gap one needs to assume an unreasonably large density of 
optical radiation. This is demonstrated in Figure Q3] where 
an agreement of IC calculations with the reported gamma- 
ray fluxes is achieved assuming an additional, although in our 
view quite unrealistic, component of optical radiation with a 
density of 140 eV cm" 3 . 

It should be noted that the problem of explanation of low 
energy gamma rays in Figure QA] is related to the large en- 



ergy of the break in the electron spectrum given by equa- 
tion ([Tol l, and correspondingly to the position of the Comp- 
ton peak which appears above 1 TeV in the spectral energy 
distribution of gamma rays. Thus, the reduction of the break 
energy down to 200 GeV could in principle solve the prob- 
lem. However, since the magnetic field in this model cannot 
exceed 15 /iG, the only way to shift the Compton peak to sub- 
TeV energies is to assume that the supernova remnant is older 
than 10 4 yr. 

An alternative solution for explanation of gamma-ray data 
within the IC models is to assume a higher density of rela- 
tivistic electrons responsible for < 1 TeV gamma rays, i.e. 
to postulate a second, low-energy electron component in the 
shell. Formally this assumption does not a priori contradict 
the observations since the additional electrons are required 
to be present below 20 TeV, i.e. in the energy range which 
is not constrained by Suzaku observations (for production of 
the lowest energy X-rays detected by Suzaku in a magnetic 
field of 15 /iG the electron energy must exceed 40 TeV). The 
results calculated under such assumption are shown in Fig- 
ure Q2] The second electron component is assumed to have 
a power-law injection spectrum with the same index as the 
first (main) component, s = 2.0, but with a high energy cutoff 
around 10 TeV. The latter is required to prevent the conflict 
with the observed gamma-ray spectrum above 1 TeV. 

3.4.2. Morphology 

A comparison between the X-ray and the TeV gamma-ray 
morphology is expected to provide information about the ac- 
celeration/emission processes. As is already shown in Fig- 
ure[8] the Suzaku XIS image revealed a significant correlation 
with that of H.E.S.S. telescopes not only in the bright struc- 
tures but also in the dim regions of the remnant. In the follow- 
ing, the two images are compared with each other on a more 
quantitative basis. 

It should be taken into account that the X-ray morphology 
can be affected by the spatial distribution of absorption col- 
umn density, A/h- In order to avoid the effect, here we disre- 
gard the energy band lower than 1 keV, and use two energy 
intervals, 2-5 keV and 5-10 keV, for comparison. For the 
A/h variation as shown in Table 3, the count rate can vary by 
5% and 0.7% in the 2-5 keV band and 5-10 keV band, re- 
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FIG. 13. — Multi-wavelength spectral energy distribution (SED) of 
RX J1713.7— 3946 with a hadronic model. The magnetic field is 200 fj,G. 
The injection index of electrons and protons is 2.0, and constant injection 
during 1000 yr is assumed. For the IC flux, the contribution of each photon 
field is shown. The dashed line, the dash-dotted line, and the dotted line are 
IC flux of CMB, infrared, and optical photons, respectively. By taking 1 kpc 
as the distance of the SNR, the total energy of electrons is W e = 3 . 1 X 1 46 erg. 



The total proton energy is W p ■■ 
ambient matter density. 



2.7 X 10 (n/l cm ) erg, where n is the 
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FIG. 14.— TheSEDofRXJ1713.7-3946 withaleptonicmodel. The in- 
jection index of electrons is 2.0. The calculation is for magnetic field strength 
of 14 fiG. The total energy of electrons is W e = 1. 4 X 10 47 erg The line styles 
for the IC spectra are the same as those in Figure [T3l 
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FIG. 15. — Same as !14l but additional optical radiation is included as seed 
photons of IC scattering to fit the H.E.S.S. spectrum. The IC spectrum with 
the thin solid line and the dotted line corr espond to the contribution of the in- 
terstellar radiation assumed in Figure [141 and the additional optical radiation, 
respectively. The thick solid line indicates the sum of the two. The energy 
density of the additional optical radiation is 14 eV cm -3 , which i s about two 
order of magnitude larger than the estimate by Porter et al. 1 2006). 



spectively. One should also take into account the difference 
of point spread functions between Suzaku XIS and H.E.S.S. 
when comparing the two images. To prevent this, we com- 
pare surface brightness for each square region with a size of 
10'. 8 x 10'. 8, which is larger than the point spread functions 




Log(E/eV) 

FIG. 16. — Same as[l4]but another electron population is added. The thin 
solid lines are radiation from each electron population and the thick solid line 
is the sum of the contributions from the two electron populations. The total 
energy of the second electron population is W e = 3.4 X 10 47 erg 

of either observatory. The regions are indicated with green 
dashed lines in Figure [TTI 

Figure [18] shows a scatter plot between the Suzaku XIS 
count rate (fk e v) and that of H.E.S.S. (Frev)- As seen in this 
plot, the X-ray count maps correlate strongly with the gamma- 
ray count map. The correlation coefficients are calculated to 
be 0.85 and 0.83, for 2-5 keV band and 5-10 keV band, re- 
spectively. It is worth noting that there are some deviations 
(X-ray intensity excesses) in the bright regions. Figure [19] 
shows a map of F^ e y-Fj e y (for 2-5 keV band) overlaid with 
the H.E.S.S. contours, in which one can see that the X-ray 
excesses are present along the NW and SW rims. 



4. DISCUSSION 



4.1. 



Cutoff in the Synchrotron Spectrum 

We conducted a series of Suzaku observations which cov- 
ers about two-thirds of the surface of SNR RX J 1 7 1 3 .7-3946. 
Through the data analysis, we successfully detected signals 
up to ~ 40 keV from each of the pointings. The HXD spec- 
tra above 10 keV are significantly steeper than those obtained 
from the XIS below 10 keV, suggesting that a spectral cutoff is 
common throughout the remnant. By combining the XIS and 
HXD spectra, we obtained a wide-band spectrum with high 
statistics, which clearly shows a cutoff around 10 keV. 

Taking advantage of the high photon statistics, we 
performed a detailed study of the cutoff shape and 
compared it with a r e cent theoretical prediction by 
IZirakashvili & Aharonianl (120071) . A sharp cutoff of the 
accelerated electron spectrum is needed to reproduce the 
cutoff shape in the synchrotron spectrum detected with 
Suzaku. The spectrum of electrons derived from Suzaku 
data is in good agreement w ith the analytical model of 
IZirakashvili & Aharonianl d2007l) . 

The cutoff energy in the spectrum of synchrotron radiation 
contains an important information about the efficiency of dif- 
fusive shock acceleration. For acceleration in the Bohm dif- 
fusion regime and when energy losses of electrons are domi- 
nated by synchrotron cooling, the cutoff energy, sp in equation 
(O is expressed as (Zirakashvili & Aharonian 2007) 



eo = 0.55 



>1 



keV, 



(11) 



v 3000 km s" 

where v s is the shock speed and 77 (> 1) is the so-called "gy- 
rofactor". The case of 77 = 1 corresponds to the "Bohm limit", 
and implies high level of turbulence 5B ~ B. The Suzaku spec- 
trum is characterized by the best-fit parameter £0 = 0.67 keV 
which gives v s = 3300/7 1 ' 2 km s" 1 . Here, we assume that 
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FIG. 17. — The regions for comparison of the X-ray and the TeV gamma- 
ray morphologies o verlaid on the XIS ima ge (2-5 keV) (upper panel) and 
the H.E.S.S. image I Aharonian et al. 2007) (lower panel). The size of each 
square region is 10'. 8 X 10'. 8. The results of the comparison are shown in 
Figure fT8l and Figure [T9l 
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FIG. 18. — The scatter plot of X-ray surface brightness from the XIS and 
the gamma rays from H.E.S.S. for each square region shown in Figure [TTl 
Each data set is scaled for direct comparison. 



FIG. 19. — The map of FkcV _ ^TcV: both of which are plotted in Figure [T8l 
The values of fjteV are those for the energy range of 2-5 keV. Th e overlaid 
contours are the H.E.S.S. image taken from Aharonian et al. ( 2007). 

the shock speed v, is uniform throughout the remnant, which 
is supported by the fact that the outer boundary of the X- 
ray morphology is nearly circular. The upper-limit of the 
shock speed v x < 4500 km s" 1 derived from the Chandra data 
( Uchiyama et al. 2007), results in ?y < 1.8. This is a strong ev- 
idence of acceleration of electrons in the regime close to the 
Bohm limit. Note that a s imilar result was obtai ned for the 
SW rim of the remnant in Uchiyama et al. (2007). Here we 
confirm this conclusion for a larger area of the remnant with 
higher statistics. 

4.2. Multi-wavelength Spectrum 

While there is little doubt in the synchrotron origin of 
broad-band X-ray emission measured by Suzaku, the X-ray 
spectrum alone does not give preference to the strength of the 
magnetic field in the region of production of synchrotron radi- 
ation. Formally, the field can be as small as 10 /iG and as large 
as 100 /iG. Meanwhile, the strength of the magnetic field has 
dramatic impact on the origin of TeV gamma-rays. The so- 
called leptonic or inverse Compton models require magnetic 
field between 10 /iG and 15 /iG. Even so, it is difficult to 
achieve, at least within a simple one-zone model, a satisfac- 
tory explanation of both X-ray and TeV gamma-ray spectral 
features, unless we invoke an extremely high diffuse radiation 
field of optical photons to enhance the IC gamma-radiation 
below 1 TeV (see Figure [15] ). A more realistic approach 
for explanation of the broad-band TeV gamma-ray spectrum 
within IC models can be realized under the assumption of ex- 
istence an additional, low-energy electron component in the 
shell (see Figure [T6l>. Even so, the most serious problem for 
IC models remains the requirement of low magnetic field in 
the gamma-ray production region, in contrast to large mag- 
netic field required to explain the fast variability of X-ray 
emission on small scales. Formally, one may assume that 
gamma-rays are mainly produced in "voids", i.e. in regions 
with very low magnetic field. This would imply quite inho- 
mogeneous distribution of the magnetic field in the shell. One 
the other hand, the observed strong X-ray and TeV correla- 
tion within the IC models can be explained only in the case of 
homogeneous distribution of magnetic field. 

The large-scale magnetic fields on parsec scales with an 
average strength larger than > 15 /iG make the IC gamma- 
ray production inefficient, and thus give preference to the so- 
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FIG. 20.— The SED of RX J1713.7-3946 from radio to X-ray with a 
synchrotron model curve calculated by assuming B = 100 fiG. 

called hadronic models of gamma-rays produced at interac- 
tions of accelerated protons with the ambient gas via pro- 
duction and decay of secondary 7r°-mesons. What concerns 
X-rays, they are produced, as in leptonic models, by syn- 
chrotron radiation of directly accelerated electrons. This is 
demonstrated in Figure Q~3] for very strong magnetic field, 
B = 200 /xG. Note that while comparing the model predic- 
tions with measurements in the radio band, one should take 
into account that the radio points shown correspond to mea- 
surements of NW rim, while the X-ray and gamma-ray points 
are for the entire remnant. If the ratio of the radio flux from 
the NW rim to that from the whole remnant is not much differ- 
ent from the corresponding ratio in X-rays, the flux from the 
whole SNR should be significantly larger. This would reduce 
the difference between the measurements and predictions. In 
any case, the radio flux can be significantly reduced assuming 
somewhat smaller magnetic field or harder electron spectrum. 
Indeed, in Figure [20] we show model calculations performed 
for a magnetic field B = 100 /iG. While the synchrotron X-ray 
flux is described perfectly as before (in Figure [13), the radio 
flux is by a factor of four lower; at 1 .4 GHz it is 34 Jy, which 
is close to the latest estimates of radio flux from the whole 
remnant based on observations with ATCA and the 30-m ra- 
dio telescopes of IRA (F. Acero et al. 2008, in preparation; G. 
Dubner 2008, private communications). 

The radio flux can be suppressed even for the ambient field 
larger than 100 /iG, provided that the electron injection spec- 
trum is harder than E~ 2 . Figure [21] demonstrates this pos- 
sibility, where we assume an electron/proton index of s = 1 .7 
which corresponds to a compression ratio of a = 5.3. Note that 
a can exceed the adiabatic upper limit of 4 as described by 
iBerezhko & Volkl(l2006l) . N ote that the value of s = 1 .7 is con- 
sistent with the conclusion of lVillante & VissanH ((2007) based 
on semi-analytical derivation of the parent proton spectrum 
from the H.E.S.S. data. In model calculations shown in Fig- 
ure [21] the spectrum of protons requires an "early" exponen- 
tial cutoff at E p o = 25.0 TeV. Note that formally the spectral 
index s = 1 .7 implies shock acceleration in non-linear regime 
which in fact predicts some deviation from pur e power-law 
distribution of accelerated par ticles (see e.g. lEllison et al.l 
I2007t IBerezhko & Volkll2006l) . This would lead to further 
reduction of the radio flux. 

The convection of low energy electrons could be another 
reason for low radio flux. Note that the escape of electrons 
through convection has strong impact only on low energy 
electrons; because of fast synchrotron cooling, the effect of 
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FIG. 21 .— The SED of RX J1713. 7-3946 with a hadronic model when the 
index of the electron/proton spectrum is 1.7. The magnetic fiel d is 2 00 fiG. 
The line styles for the IC spectra are the same as those in Figure [73l the total 
energy of electrons is W e = 6.0 X 10 45 erg. The total proton energy is W p = 
1.6 X 10 50 (n/l cnr 3 r' erg. 

escape is negligible for multi-TeV electrons. A significant 
quantity of low-energy electrons can escape from the shell of 
the SNR before emitting radio photons. Therefore, the ra- 
dio flux can be reduced while the X-ray flux will remain un- 
changed. 

For calculations shown in Figure Qj] the to- 
tal energy of electrons is estimated as W e = 3.1 x 
10 46 (d/1 kpc) 2 erg, and the energy for protons as 
W p = 2.7 x 10 50 (n/l cnT 3 rV/l kpc) 2 erg. In the case of 
harder energy spectra with power-law index s = 1.7 corre- 
sponding to Figure[2T] one has W e = 6.0 x 10 45 (d/l kpc) 2 erg, 
and W p = 1.6 x 10 50 (n/l cm~ 3 ) _1 (c//l kpc) 2 erg, respec- 
tively. The proton/electron ratio in either case is very small, 
K ep < 10" 4 (n/l cm" 3 ). This value is significantly smaller 
than that for directly observed local cosmic rays (K ep ~ 0.01), 
unless a la rge ambient matt e r dens ity of n ~ 100 c m" 3 is 
assumed. iKatz & Waxmanl d2008l) and iButt et all d2008l) 
argued that the hadronic scenario for this SNR has difficulties 
because the K ep value should be consistent with local cosmic 
rays and other SNRs. However, the K ep value of one SNR 
at a fixed age does not necessarily need to agree with the 
local cosmic ray value. The low-energy electrons are likely 
produced in later stages of SNR evolution, when the K ep 
value can be different from the present value. A comparison 
with other SNRs should be performed with care, as well. 

Cutoff energy in the gamma-ray spectrum should give an 
important hint whether SNRs are sources of cosmic rays 
below the knee, if gamma ray s observed by H.E.S.S. have 
hadronic origins. Plaga (2008) argued that the cutoff energy 
of the HE. S .S. spectrum of RX J 1 7 1 3-3946 is around 1 8 TeV, 
which can be translated to proton energy more than 10 times 
below the energy of the knee. Indeed our multi-wavelength 
study requires a cutoff in the proton spectrum around 100 TeV 
or even less for hard acceleration spectra of protons. However, 
one should take into account that even in the case of effective 
acceleration the highest energy protons beyond 100 TeV es- 
cape the source in a quite short time scales, and hence do not 
contribute to the gamma-ray production at the present ep och 
dPtuskin & Zirakashvilill2 005; Gab ici & Aharonianll2007[) . 

A unique feature of RX J 17 13.7-3946 is the l a ck of ther- 
mal X-ray emiss ion. Recently, IKatz & Waxmanl (120081) and 
IButt et ai] d2008l) interpreted this fact as an argument against 
the hadronic model for TeV gamma rays. Generally it is true 
that plasma in young supernova remnants is heated to high 
temperatures observed via thermal X-ray emission of hot elec- 
trons. However, one should take into account that we deal 
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with a unique object, and the lack of thermal X-ray emis- 
sion cannot a priori be invoked as an argument against the 
hadronic origin of the observed TeV gamma rays. 

It is important to note that in SNR shocks the formation 
of high plasma temperatures with kT, ■ = 3/16 m.{V* is relevant 
only to protons (ions), and that a high ion temperature does 
not automatically (from first principles) mean a high electron 
temperature. In fact, the only known heating process of ther- 
mal electrons is Coulomb collisions between electrons and 
protons (ions), which, however, has too long time scale to es- 
tablish electron-proton equipartition. On the other hand, we 
do know from X-ray observations that the electrons in young 
SNRs are heated to keV temperatures. This can be explained 
by assuming that a hypothetical mechanism, most likely re- 
lated to the energy exchange through excited plasma waves, is 
responsible for effective electron heating in SNRs. As long as 
the nature of this mechanism in collisionless shocks remains 
unknown, one cannot predict, even qualitatively, the specifics 
of its operation on a source by source basis. 

We indeed deal with two interesting facts. First, many 
young SNRs, like Tycho and Cassiopeia A, with intense ther- 
mal X-ray emission and intense nonthermal radio emission 
emit little (or do not emit at all) TeV gamma rays. On the 
other hand, RX J1713. 7-3946, with lack of (or rather very 
low) thermal X-ray emission and with relatively weak non- 
thermal radio emission, is a source of powerful TeV radiation. 
These two facts can be treated as a hint for low efficiency of 
establishing equipartition in the thermal plasma in very ef- 
fective TeV particle accelerators like RX J1713. 7-3946. It 
is interesting to note in this regard, that such a tendency is 
found also for an other effective TeV accelerator — SN 1006 
dVink et al.l2003l) . Whether the reduction of the exchange rate 
between different particle species in thermal plasma has a link 
to the partic le accelerat i on in high Mach number shocks, as 
proposed by IVink et al.l (12003), is a very interesting question 
to be explored in future deep theoretical and phenomenolog- 
ical studies. In this regard, RX J1713. 7-3946 can serve as a 
k ey "template" s ource for such studies. 

Hug hes et al.l (2000) also discussed low electron tempera- 
ture based on Chandra observations of a young SNR in the 
Small Magellanic Could, IE 0102.2-7219. They measured a 
blast-wave velocity of ~ 6000 km s" 1 from the expansion rate 
and predicted the electron temperature of k T e > 2.5 keV by 
considering Coulomb heating. However, the electron temper- 
ature derived from their spectral analysis is 0.5 keV, which is 
far below the prediction. According to their discussion, not 
only electron heating but also ion heating is suppressed and 
substantial fraction of energy may be going into cosmic -ray 
production due to the non-linear effects in the shock. 

The nonlinear shock acceleration in this object can con- 
vert a significant, up to / ~ 0.5 fraction of the kinetic energy 
of explosion into relativistic particles. Correspondingly the 
fraction of available energy which goes to the heating of the 
ambient plasma will be reduced 1 — / ~ 0.5. Yet, conserva- 
tive estimates show that plasma in RX J1713. 7-3946 can be 
heated to quite high temperatures even the heating of elec- 
trons and protons proceeds only through the Coulomb ex- 
chang e. This question recently has studied by Ellis on et alJ 
(2007) for a standard SNR of age ?snr = 500 yr and energy 
£sn = 10 51 erg. In particular, it has been shown that in the 
case of effective diffusive shock acceleration and the plasma 
density n = 0.1 cm -3 , the ratio of synchrotron luminosity to 
thermal (bremsstrahlung) luminosity can be as large as 100. 
This implies that in the case of RX J 17 13.7-3946, from which 



thermal X-ray emission is not observed, the plasma density 
cannot significantly exceed 0.1 cm" 3 (the luminosity of ther- 
mal bremsstrahlung is proportional to n 2 ). The low den- 
sity of order of 0.1 cm" 3 reduces the parameter space for 
hadronic models but does not exclude it. Indeed, as men- 
tioned above, the TeV gamma-ray flux of RX J1713. 7-3946 
can be explained by interactions of protons if the parameter 
A = (Wp/10 50 erg)(n/l cm" 3 )(t//l kpc)" 2 exceeds 1.5 to 3, de- 
pending on the spectrum of protons. Assuming that more than 
30% of the explosion energy of this SNR is released in accel- 
erated protons, and that the plasma in the gamma-ray pro- 
duction region is compressed by a factor of few, we find that 
the located of the source at a distance of about 1 kpc would 
marginally support the hadronic model. While closer location 
of the source would make the model requirements quite viable 
and flexible, the location of the source beyond 1 kpc hardly 
can be accommodated within a standard shock acceleration 
scenario. If the SNR is closer, the distance of d > 0.5 kpc 
seems reasonable considering the upper limit on the sh ock 
speed of 4500(<i/l kpc) km s" 1 by lUchiyama et all (120071) . 

Finally we shou ld mention the model suggested by 
Malk ovet alJ (|2005) which can naturally explain both the low 
synchrotron flux at radio frequencies and lack of thermal X- 
ray emission of RX J1713. 7-3946. The standard scenarios 
of gamma-ray production in SNRs assume that radiation is 
produced in downstream where the density of both relativis- 
tic particles and thermal plasma is higher than in upstream. 
However, in the cases when the shock is expanding into a low- 
density wind bubble and approaching cold dense material, e.g. 
swept-up shell or surrounding molecular clouds, the gamma- 
radiation is contributed predominantly from upstream. While 
the energy distribution of accelerated particles downstream is 
coordinate-independent in both linear and nonlinear regimes, 
the particle distribution upstream is coordinate-dependent. 
Because of energy-dependent diffusion coefficient, the high- 
energy particles diffuse ahead of low-energy particles, thus a 
dense material adjacent upstream will "see" relativistic parti- 
cles (protons and electrons) with low-energy exponential cut- 
off, E m \ n , which depends on the location of the dense regions. 
This implies that the effective production of TeV gamma-rays 
(from p-p interactions) and X-rays (from synchrotron radia- 
tion of TeV electrons) will be not accompanied by low energy 
(GeV) gamma-rays and synchrotron radio emission. Obvi- 
ously, this model is not constraint by lack of thermal emis- 
sion. 

4.3. Morphology 

In addition to the spectral information, the comparison of 
X-ray and TeV gamma-ray images presented in Figure [18] 
helps us to draw the physical picture of RX J1713.7-3946. 
Let us first discuss the tight correlation observed in most parts 
of the remnant. Within the hypothesis of hadronic origin of 
gamma-rays, the gamma-ray flux is proportional to the num- 
ber densities of the ambient matter and relativistic protons, 
while the X-ray flux is proportional to number density of elec- 
trons. If the matter distribution significantly varies throughout 
the SNR, we need fine parameter tuning among the matter dis- 
tribution, the electron injection rate, and the proton injection 
rate, in order to produce the tight correlation. Therefore, more 
natural explanation is that the matter density is uniform and 
the injection rate of the electrons and that of the protons are 
proportional to each other. 

The X-ray flux excess along the NW and SW rims provides 
a unique probe of recent acceleration activity. Let us consider 
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FIG. 22. — The same plot as Figure [T8l but results obtained by assuming the 
"toy" model described in the text. 

a "toy" model and compare its predictions with the observa- 
tional results. In the toy model, the injection rate of electrons 
and protons keeps constant but increases by a factor of 1 .5 in 
the last 10 yr only at the NW and SW rims. What is important 
here is the difference of cooling time between electrons and 
protons. The synchrotron cooling time of an electron emitting 
synchrotron photons with energy of e is given as equation ©, 
while the cooling time of protons due to p-p interactions is 
expressed as 



f D „ = 5.3 x 10 7 



1 cm 



yr 



(12) 



and is almost energy-independent. For the magnetic field 
B = 100 fiG, the cooling times of electrons emitting 2 keV 
and 5 keV X-rays are 12 yr and 7.6 yr, respectively. For any 
reasonable density of ambient gas, the cooling time of pro- 
tons is much longer than the lifetime of this SNR. While the 
synchrotron X-rays we observe at present are emitted by elec- 
trons accelerated during the last ~ 10 yr. the flux of 7r°-decay 
gamma rays is provided by protons accelerated throughout the 
lifetime of the SNR. Figure [22] shows the scatter plot of F^y 
and Fxev expected from the toy model, which shows a simi- 
lar distribution to the observational results in Figure [18] The 
recent active acceleration increases the X-ray flux while keep- 
ing the gamma-ray flux almost unchanged. 

5. SUMMARY 

We observed SNR RX J1713.7-3946 with the Suzaku ob- 
servatory. Hard X-rays up to ~ 40 keV are detected from 
each of the 11 pointings. The hard X-ray morphology es- 
timated by the HXD PIN data is generally consistent with 
an extrapolation of that in the energy region below 10 keV. 
When the HXD spectra are fitted with a power law, the pho- 
ton indices are larger than those obtained from the XIS data. 
The difference of photon indices between the XIS and the 
HXD varies from region to region. Although this may sug- 
gest a variation of spectral shape or cutoff energy, the FoV 
of the HXD PIN is not small enough to allow detailed imag- 
ing spectroscopy. Such studies will become possible with up- 
coming missions with hard X-ray mirrors like NeXT, NuS- 
TAR, or Simbol-X. Moreover these missions will extend the 



effective studies up to 100 keV. If the TeV gamma-ray emis- 
sion of RX J1713. 7-3946 is of hadronic origin, this should 
allow detection and study of synchrotron X-ray emission of 
electrons produced in proto n-proton interact ions via decays 
of secondary charged pions (lAharon ian 2004). In this regard, 
the above mentioned missions can provide very effective tools 
for deep (at the level as low as 10~ 13 erg cm s" 1 ) probes of 
hadronic processes in SNRs in the > 100 TeV energy regime 
with unprecedented (subarcmin) angular resolution. 

Using the XIS and HXD data, we obtained a synchrotron 
spectrum in the energy range of two decades (0.4-40 keV), 
which means we can probe the parent electron distribution 
in the energy range of one decade regardless of the mag- 
netic field strength. The wide-band coverage enables us to 
see a clear high-energy cutoff in the synchrotron spectrum, 
and for the first time allows detailed studies of the cutoff 
shape and derivation of the spectrum of the parent electrons. 
The spectral shape in the cutoff region was quantitatively 
evaluated, which revealed that the cutoff shape is compati- 
ble wi th a theoretical prediction by Zirakash vili & Aharonianl 
(2007). Based on their model, the cutoff energy is obtained as 
£o = 0.67 ± 0.02 keV. This result, t ogether with upper l imit of 
the shock velocity from Chandra ( Uchiyam a et alJl2007l) . in- 
dicates that acceleration efficiency in this remnant approaches 
the Bohm limit. 

We modeled the multi-wavelength spectrum of 
RX J 17 13.7-3946 within both the leptonic and hadronic 
gamma-ray production scenarios. The hadronic model 
nicely fits the data while it seems difficult to explain all 
observational data with a simple leptonic model. The major 
problem with leptonic models is related to the requirement 
of low magnetic field, B < 15 pG which does not agree with 
the recent discovery of variability of local regions of the shell 
on year timescales. The hadronic models recently have be 
criticized based on the lack of thermal emission and low radio 
flux. However, these fact can be accommodated within a 
standard shock acceleration model in a young SNR with mag- 
netic field as large as 100 pG, plasma density n ~ 0.1 cm , 
and distance to the source d < 1 kpc, assuming that more 
than 30% of energy of explosion is released in relativistic 
protons with spectrum as hard as E~ 2 and exponential cutoff 
around 100 TeV. 

Finally, besides the general strong correlation between X- 
ray and TeV gamma-ray emission, we found an excess X-ray 
emission compared to gamma-ray emission in the brightest re- 
gions of the remnant. This excess can be explained by recent 
activity accompanied by effective acceleration of electrons in 
localised reagions of the shell. 



The authors would like to thank all the members of the 
Suzaku Science Working Group for their help in the space- 
craft operation, instrumental calibration, and data process- 
ing. The authors thank Una Hwang for carefully reading the 
manuscript and Stefan Funk for assisting with the comparison 
of the morphologies with the H.E.S.S. data and for his helpful 
comments on the manuscript. The authors also thank Misha 
Malkov and Vladimir Zirakashvili for discussions related to 
different aspects of the diffusive shock acceleration theory. T 
Tanaka and A. Bamba are supported by research fellowships 
of the Japan Society for the Promotion of Science for Young 
Scientists. 



16 



Tanaka et al. 



REFERENCES 



Aharonian, F. A., et al. 2001, A&A, 370, 1 12 

Aharonian, E A. 2004, Very High Energy Cosmic Gamma Radiation: A 
Crucial Window on the Extreme Universe (World Scientific) 

Aharonian, F. A., et al. 2004, Nature, 432, 75 

Aharonian, F. A., et al. 2006, A&A, 449, 223 

Aharonian, F. A., et al. 2007, A&A, 464, 235 

Berezhko, E. G., & Volk, H. J. 2006, A&A, 451, 981 

Bird, A. J., et al. 2006, ApJ, 636, 765 

Blandford, R., & Eichler, D. 1987, Phys. Rep., 154, 1 

Butt, Y, Porter, T., Katz, B., & Waxman, E. 2008, MNRAS, in press 

Cassam-Chenai, G., Decourchelle, A., Ballet, J., Sauvageot, J.-L., Dubner, 
G., & Giacani, E. 2004, A&A, 427, 199 

Ellison, D.C., Patnaude, D.J., Slane, P., Blasi, P., Gabici, S. 2007, ApJ, 661, 
879 

Enomoto, R., et al. 2002, Nature, 416, 823 

Fukui, Y., et al. 2003, PASJ, 55, L61 

Gabici, S., & Aharonian, F. A. 2007, ApJ, 665, L131 

Gruber, D. E. Matteson, J. L.; Peterson, L. E., & Jung, G. V. 1999, ApJ, 520, 
124 

Hartman, R. C, et al. 1999, ApJS, 123, 79 

Hiraga, J. S., Uchiyama, Y, Takahashi, T., & Aharonian, F. A. 2005, A&A, 
431,953 

Hughes, J. P., Rakowski, C. E., & Decourchelle, A. 2000, ApJ, 543, L61 
Ishida, M., et al. 2006, Suzaku Memo, JX-ISAS-SUZAKU-MEMO-2006-40 
Ishisaki, Y, et al. 2007, PASJ, 59, SI 13 

Katz, B., & Waxman, E. 2008, Journal of Cosmology and Astroparticle 
physics, in press 

Kelner, S. R., Aharonian, F. A., & Bugayov, V. V. 2006, Phys. Rev. D, 74, 
034018 

Kokubun, M., et al. 2007, PASJ, 59, S53 

Koyama, K., Petre, R., Gotthelf, E. V., Hwang, U., Matsuura, M., Ozaki, M., 

& Holt, S. S. 1996, Nature, 378, 255 
Koyama, K., Kinugasa, K., Matsuzaki, K, Nishiuchi, M., Sugizaki, M., Torii, 

K, Yamauchi, S., & Aschenbach, B. 1997, PASJ, 49, L7 
Koyama, K., et al. 2007, PASJ, 59, S23 

Lazendic, J. S., Slane, P. O., Gaensler, B. M., Plucinsky, P. P., Hughes, J. P., 
Galloway, D. K, & Crawford, F. 2003, ApJ, 593, L23 



Lazendic, J. S., Slane, P. O., Gaensler, B. M., Reynolds, S. P., Plucinsky, P. 

P., & Hughes, J. P. 2004, ApJ, 602, L201 
Malkov, M.A., Diamond, PH., Sagdeev, R.Z. 2005, ApJ, 624, L37 
Malkov, M. A., and Drury, L. O'C. 2001, Rep. Prog. Phys., 64, 429 
Mitsuda, K., et al. 2007, PASJ, 59, S 1 

Mizuno, T., et al. 2006, Suzaku Memo, JX-ISAS-SUZAKU-MEMO-2006-42 
Moriguchi, Y, Tamura, K., Tawara, Y, Sasago, H., Yamaoka, K., Onishi, T., 

& Fukui, Y. 2005, ApJ, 631, 947 
Muraishi, H, et al. 2000, A&A, 354, L57 

Pfeffermann, E., & Aschenbach, B. 1996, in Rontgenstrahlung from the 
Universe, ed. H. U. Zimmermann, J. Trumper, & H. Yorke, MPE Rep., 
263, 267 

Plaga, R. 2008, Nature, 453, 48 

Porter, T. A., Moskalenko, I. V., & Strong, A. W. 2006, ApJ, 648, L29 
Ptuskin, V. S., & Zirakashvili, V. N. 2005, A&A, 429, 755 
Serlemitsos P. J., et al. 2007, PASJ, 59, S9 

Slane, P., Galensler, B. M., Dame, T. M., Hughes, J. P., Pluncinsky, P. P., & 

Green, A., 1999, ApJ. 525, 357 
Slane, P., Hughes, J. P., Edgar, R. J., Plucinsky, P. P., Miyata, E., Tsunemi, 

H, & Aschenbach, B. 2001, ApJ, 548, 814 
Takahashi, T., et al. 2007, PASJ, 59, S35 
Takahashi, T, et al. 2008, PASJ, 60, S131 (Paper I) 
Terada, Y, et al. 2005, IEEE Trans. Nucl. Sci., 52, 902 
Uchiyama, Y, Aharonian, F. A., & Takahashi, T. 2003, A&A, 400, 567 
Uchiyama, Y, Aharonian, F. A., Takahashi, T, Hiraga, J. S., Moriguchi, Y, 

& Fukui, Y. 2005, High Energy Gamma-Ray Astronomy, 745, 305 
Uchiyama, Y, Aharonian, F. A., Tanaka, T., Takahashi, T, & Maeda, Y. 2007, 

Nature, 449, 576 
Villante, F. L., & Vissani, F. 2007, Phys. Rev. D, 76, 125019 
Vink, J., Laming, J. M., Gu, M. F, Rasmussen, A., & Kaastra, J. S. 2003, 

ApJ, 587, L31 

Wang, Z. R., Qu, Q.-Y, & Chen, Y. 1997 A&A, 318, L59 
Watanabe, S., et al. 2007, Suzaku Memo, JX-ISAS-SUZAKU-MEMO-2007- 
01 

Zirakashvili, V. N., & Aharonian, F. 2007 A&A, 465, 695 



